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1. 


I.  INTKOimmON 

The  power  spent  ra  of  I Inoar  1<TI  sonar  data  for  two  different 
submarine  aspects  have  been  computed.  These  results  are  eomtvuvd 
with  each  other  and  with  t tie  power  spectra  of  the  associated  rover 
bo  rat  l on  dut  a . 


In  conjunct  ion  w i t ti  t tie  t'lpnal  I'lyvslcs  Division  of  t lie  .Pe  fenne 
Keneaivh  L'lhorat  ory  ( PKt  a different  technique  t’or  d i a pi  ay  i nr. 
analyzed  sonar  data  lias  been  Investigated.  Tills  display  consists  of 
quadrature  sampled  data  and  eorivlations  t tint  are  plotted  with  t he 
eomputer  as  project  ions  in  t tiree  dimensions. 

A meet  Ins;  was  hold  between  the  Texas  Instruments,  Inc.  (Ti), 
Kqulimient  hraneti  and  t tie  i'Kl.  Computer  Science  Division.  The  sul'.leet 
of  tlw'  meeting  was  nonce  rued  with  application  of  Tl’s  met  hod  ot' 
curve  fitting  in  inviie  types  of  data  analysis. 
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II.  POWER  SPECTRA 

((?)  In  order  to  gain  a better  intuitional  feel  for  the  analysis  of 
actual  sonar  data,  several  power  spectra  were  computed  for  two 
different  submarine  aspects.  These  spectra  were  computed  by  squaring 
the  modulus  of  a discrete  finite  Fourier  transformation.  The  trans- 
mitted signal  was  a linear  FM  upsweep  with  a bandwidth  W = 625  Hz, 
pulse  length  T = 0.1024  sec,  TW  = 64,  and  the  return  data  was  sampled 
at  four  times  the  carrier  frequency,  Fq  = 5000  Hz.  Only  the  frequency 
band  of  interest  was  plotted,  and  each  scale  was  normalized. 

Figures  1 and  2 show  the  power  spectra  for  a beam  aspect  submarine. 

Of  particular  interest  is  the  low  frequency  ripple  in  the  power 
spectra.  This  is  the  information  that  Cepstrum  and  other  correlation 
techniques  use  to  decompose  the  data  into  highlights.  As  an 
interesting  calculation,  the  distance  between  the  large  peaks  in 
Fig.  1 is  about  122  Hz.  This  information  can  be  related  to  a physical 
distance  as 


The  beam  of  this  particular  submarine  is  approximately  27  ft.  Even 
though  this  was  a rough  calculation,  the  numbers  are  of  the  same 
magnitude.  The  purpose  of  the  study  was  to  demonstrate  why  the  power 
spectra  of  the  stern  aspects  shown  in  Figs.  3 and  4 are  so  much 
different  than  the  spectra  of  the  beam  aspects.  This  is  to  be 
expected  because  the  highlights  of  a stern  aspect  would  be  separated 
more  in  time  than  the  highlights  for  a beam  aspect.  Since  the 
frequency  domain  is  inversely  proportional  to  the  time  domain,  the 
ripple  in  the  spectra  for  a stern  would  be  much  more  closely  spaced 
than  the  ripple  in  the  spectra  for  a beam  aspect.  The  characteristics 
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POWER  SPECTRA  OF  STERN  ASPECT  SUBMARINE  DATA  (U) 
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(C)  of  the  power  spectra  for  the  reverberation  shown  in  Figs . 5 and  6 
•ire  not  much  different  than  those  for  the  stern  aspect.  It  is  felt 
that  the  peaks  in  the  reverberation  spectra  will  be  randomly  spaced 
and  will  not  transform  back  to  give  apparent  highlights,  but  this 
remains  to  be  seen. 
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III.  THREE-DIMENSIONAL  PLOTS 


(C)  For  active  sonar  a potentially  powerful  tool  in  signal  processing 
is  the  detection  of  phase  in  returning  highlights  with  respect  to  the 
transmitted  signal.  This  phase  information  can  be  used  in  addition 
to  locating  highlights  in  time  and  measuring  their  amplitudes.  As 
reported  in  Quarterly  Progress  Report  No.  5 on  Contract  NObsr-95l8l  (u) 
it  was  demonstrated  on  artificial  data  that  Medgyessy  analysis  (impulse 
response)  would  compute  the  phase,  location  in  time,  and  amplitudes 
of  highlights . However,  with  real  data  the  present  method  of  com- 
puting a Medgyessy  analysis  gives  poor  results.  So  far  the  results 
have  been  very  fuzzy,  looking  much  like  white  noise.  This  is 
attributed  to  the  fact  that  narrow  band  data  are  treated  as  wide  band 
data.  The  assumption  is  that  noise  outside  the  band  is  masking  the 
expected  result.  It  is  felt  that  a smoothed  Medgyessy  analysis  can 
be  computed,  but  this  has  not  been  tried. 

(U  -FOUG^  Another  way  of  approaching  the  problem  of  phase  detection  is  to 

* 

process  the  data  in  quadrature.  With  this  type  of  processing,  it 
is  suggested  that  the  results  be  displayed  in  three  dimensions,  x(t), 
y(t),  and  time  t.  As  an  example,  consider  the  two  plots  in  Fig.  7. 
These  plots  were  computed  from  artificial  data  where  a linear  PM 
upsweep  was  used  as  a transmitted  signal.  The  difference  between  the 
two  plots  is  the  TW  product.  The  transmit  signal  was  crosscorrelated 


0.  D.  Grace  and  S.  P.  Pitt,  "Quadrature  Sampling  of  High  Frequency 
Waveforms,"  to  be  published  in  correspondence  section  of  The  Journal 
of  the  Acoustical  Society  of  America. 


**. 


Narrow  band  signal  s(t)  = x(t)  cos u>  t + y(t)sin  u^t;  x and  y are 
the  quadrature  components  and  is  the  carrier  frequency. 


<i  \ r w,  m**  ;t«  *• 
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with  four  returning  highlight  a,  when*  the  third  highlight  waa 
ItV  deg  out  of  phase  with  the  transmitted  signal.  These  curves 
demonstrate  that  the  quadrature  components  of  a crosseorrelatiou 
plotted  In  three  dimensions  will  Indicate  phase,  location  In  time, 
and  amplitudes  of  highlights. 

(0)  Quadrature  sampled  c rosscorrelat ions  were  computed  for  some  bow 
aspect  data.  The  results  are  shown  In  Fig.  8.  It  appears  that  the 
second  and  third  returns  are  upside  down  (shifted  180  deg)  with 
respect  to  the  first  return,  .lust  how  this  happened  is  not  certain. 

It  Is  known  that,  the  target  was  snorkel  ing  at  a range  of  yd, 

and  that  the  sea  surface  was  glassy  calm  in  this  ('articular  sea  test. 

A pure  conjecture  is  that  the  two  returns  are  acoustic  mirror  Images 
of  the  target  from  the  surface,  which  would  he  shifted  180  deg. 

(C)  Another  how  aspect  Is  shown  In  Fig.  *>,  and  consists  of  the 
quadrature  sampled  return  and  the  crosscorrelat Ion.  These  were 
sampled  at  a lower  rate,  and  there  Is  noticeable  looping  In  the 
plots.  The  looping  In  the  oorre.lat  Ion  looks  a little  suspicious, 
and  one  Is  always  plagued  by  the  thought  that  some  obscure  progrtunmlng 
error  has  crept  Into  the  results.  The  point  to  be  made  by  Fig.  *'  is 
that  tiie  quadrature  sampled  return  may  produce  a characteristic 
picture,  but  the  results  are  difficult  to  interpret  . However,  t lie 
quadrature  sampled  crosseorrelatiou  Is  "easy"  to  interpret.  In 
particular,  each  main  "blip"  represents  a hlglQ  tglvt  ; its  angle  with 
respect  to  the  y-axls  represents  Its  phase,  and  its  position  on  the 
t-axls  represents  its  loeat Ion  In  time.  It  is  clear  that  there  is 
phase  distortion  of  the  highlights  in  the  real  data.  That  is,  some 
of  the  "blips"  are  neither  straight,  up  nor  straight  down  (0  deg  or 
180  deg).  Tills  Is  difficult  to  see  on  a two-dimensional  surface; 
some  stereoscopic  views  (not  shown)  were  made  of  the  returns  and  of 
the  correlations.  These  wen*  made  by  producing  a plot  similar  to 
Fig.  9 in  red  color,  rotating  the  time-axis  *'  deg  about  live  y-axls, 

1 1 
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three  quadrature  sampled  crosscorrelations 
for  bow  aspect  submarine  (U) 
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QUADRATURE  SAMPLED  RETURN 


QUADRATURE  SAMPLING  (U) 
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and  plotting  a second  projection  In  green  color.  When  viewed 
through  color  filters  (red  for  one  eye  and  green  for  the  other), 
the  result  appears  as  a single  black  line  projected  above  the  paper 
background.  With  a little  practice,  the  data  can  be  seen  spiraling 
around  the  time-axis,  and  In  theory  at  least,  one  can  see  phase 
distortion,  amplitude,  and  location  la  time,  all  in  one  picture. 

( I f It  Is  felt  that  displaying  sonar  ret  urns  In  three  dimensions 
would  be  of  value,  then  In  practice  It  would  be  an  easy  thing  to  do 
For  exiunple,  a color  television  could  be  used  to  project  the  two 
stereoscopic  views  In  different  colors  and  then  the  results  viewed 
through  different  colored  lenses.  This  Is  well  within  the  capa- 
bilities of  present  technology. 
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IV.  LAG  WINDOWS  TOR  CEPSTRUM 


In  vluurterly  Progress  Report  No.  o on  Contract  N0bsr-O'jl8l  (U), 
it  was  stated  that  the  baste  shape  of  the  lag  window  applied  to  the 
autocorrelat ton  dominates  the  shape  of  the  Cepstrum.  A good  solution 
to  this  problem  Is  found  to  be  a clipped  cosine  bell. 


1 To  see  the  significance  of  this.  Fig.  10  shows  a lagged  CYpst 


omputed  with  three  different  lag  windows.  Window  11ANN  111  is 


1/0(1  t cos  .it  -T  ) 
x m 


0 < t v T 
— — in 


whore  T Is  the  greatest  Ihk  time.  The  result  is  certainly  smoothed, 
m 


but  a VI  of  the  Interesting  highlights  were  wiped  out  . The  HANN  it 


i,/o(i  * cos  t At  'i 


0 v t o'  t 


The  result  Is  a stretched  HANN  111;  but  the  autocorrelat. Ion  was 


truncated  at  T . and  this  Introduces  a spike  into  the  Cepstrum  which 
m 


HANN  '10  does  not  smooth  enough.  The  dipped  cosine  bell  is 


0 s.  t v 0.0  T 
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which  means  that  the  last  10  percent  of  the  autocorrelat- ton  Is 
smoothed,  and  the  remaining  'V*  percent  is  unbiased.  Tilts  gives  a 
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LAGGED  CEP STRUM 


UANN  111 


CLIPPED  COSINE  BELL 


0.0  TIME  - sec  0.05 


FIGURE  10 

EXAMPLES  OF  THREE  WINDOWS  USED  IN  LAGGED  CEPSTRUM  (U) 

SEQUENCE:  3 

TAPE:  147 
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good  result  in  t ho  Oepst  nan.  The  highlight  s are  unbiased.  and  t he 

spike  it  the  point  of  t runcat ton  T is  smoothed  to  zero. 

m 


Using  t ho  clipped  cosine  window , some  t lino  varying  Oepst rums 
wore  computed  for  bow  aspect  data,  'l'ho  envelopes  of  the  Oepst  nan 
wore  computed  and  are  plotted  in  Figs.  11.  1.',  I'.  lU,  and  l1'.  The 
original  data  records  of  each  return  consisted  of  <\VO  points 
ssunpled  at  four  times  the  carrier  froquonev  V ■ kll:.  The  trans- 
mitted signal  was  a linear  W,  upsweep  with  a bandwidth  W ' and 

TW  V.  v^n  each  return,  five  overl  apping  blocks  of  data  were  taken 
and  a lagged  Oepst  nun  was  comput  od  for  each,  block.  Fach  block 
consisted  of  .'uOo  points,  and  1‘aV  lags  wore  t akon  for  each  Oepst rum. 
The  snmpl  Ing  epoch  of  each  block  was  shifted  b\  W'  points,  starting 
at  point  number  one.  The  time  between  points  is  ITT  1 o'  AWV  l!.-. 
Pherefore.  aV  points  represent  kV  .'AW'  He  0.00''  sec. 

L'ho  most  significant  feature  of  this  analysis  is  tv'  show  that 
the  Oepst rum  results  are  almost  independent  of  the  sampling  epoch. 

The  target  rosolut i on  in  all  five  returns  started  around  point 
number  l uV.  If  the  sampling  epoch  starts  within  aV  points  or 
OoV'  si'C  of  the  target  resolution,  then  the  Oepst  rum  results  are 
about  the  same.  This  is  a g.Oi\t  feature  of  Oepst  rum  .and  gives  it  an 
advantage  over  other  types  of  data  processing  that  depend  cn  a 
sampling  epoch  coinciding  with  target  resolution  for  pattern 
reeogni t Ion . 
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FIGURE  11 

TIME- VARY  INC  LAGGED  CEPSTRUM  ENVELOPE 
FOR  BOW  ASPECT  SNORKELING  SUBMARINE  (U) 

SEQUENCE:  3 

TAPE:  197 
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W - 625  Hz 

F ■ 5 kHz 

o 

LINEAR  FM  UPSWEEP 

DATA  LENGTH:  2596  POINTS 
max  LAG:  1500 
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FIGURE  12 

TIME-VARYING  LAGGED  CEPSTRUM  ENVELOPE 

FOR  BOW  ASPECT  SNORKELING  SUBMARINE  U1'* 

SEQUENCE:  5 

TAPE:  197 

TW  - 32 

W - 625  Hz 

F - 5 kHz 

o 

LINEAR  FM  UPSWEEP 

DATA  LENGTH:  2596  POINTS 
max  LAG:  1500 
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TIME- VARYING  LAGGED  CEPSTRUM  ENVELOPE 

FOR  BOW  ASPECT  SNORKEL  I NG  SUBMARINE  (U> 

SEQUENCE:  7 
TAPE:  197 

rw  - j2 
w - 625  Hz 

F - 5 kHz 

o 

LINEAR  FM  UPSWEEP 
DATA  LENGTH:  25»J6  POINTS 
max  LAG:  1500 
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TIME-VARYING  LAGGED  CEPSTRUM  ENVELOPE 

FOR  BOW  ASPECT  SNORKELING  SUBMARINE  (U) 

SEQUENCE:  9 

TAPE:  197 

TW  - 32 

W - 625  Hz 
F - 5 kHz 


LINEAR  FM  UPSWEEP 
DATA  LENGTH:  2596  POINTS 
max  LAG:  1500 
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TIME-VARYING  LAGGED  CEPSTRUM  ENVELOPE 

FOR  BOW  ASPECT  SNORKELING  SUBMARINE  (U) 

SEQUENCE:  II 

TAPE:  197 


LINEAR  FM  UPSWEEP 

DATA  LENGTH:  2596  POINTS 
max  LAG:  1500 
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V.  MEETING  Win:  TEXAS  INSTRUMENTS,  INC. 

(U-J£Q110)  On  29  May  1968  a technical  meeting  was  held  between  the  Texas 

Instruments'  Equipment  Branch  and  the  DRL  Computer  Science  Division. 
TI  was  represented  by  Messrs.  MacDonald  and  Rice  and  Dr.  Spitznoggle. 
DRL  was  represented  by  Messrs.  Ellis  and  Shooter  and  Miss  Webb. 

(U-SOW)  The  purpose  of  the  meeting  was  to  discuss  applications  of 

curve  fitting  techniques  that  have  been  implemented  by  TI . As  an 
introduction  to  DRL,  Mr.  Rice  agreed  to  make  an  analysis  of  some 
unequally  spaced  time-series  data  supplied  by  Miss  Webb.  The  data 
represent  the  light  modulation  of  stars  and  is  of  no  interest  to 
the  Navy  except  as  a demonstration  of  a type  of  analysis. 

(U-POVQ*)  Since  DRL  is  interested  in  pattern  recognition,  it  is  suggested 
by  TI  that  if  an  equation  were  computed  describing  a spectrum  or  a 
correlation  (pattern),  then  it  might  be  easier  to  detect  differences 
between  equations  than,  for  example,  correlations. 

(U  DOUO“)  These  ideas  will  be  pursued  further,  and  the  results  of 
Mr.  Rice's  analysis  will  be  reported  later. 
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San  Diego  Division 
• VI  Catalina  B.lvd. 

San  Diego,  California  (VP '2 
At  in:  Code  Do06 

Coiiunanding  Officer  and  Director 

U.  S.  Nava  Underwater  Sound  Laboratory 

Fort  Trumbul I 

New  London,  Connect  lent  iWO 
Attn:  Code  224:' 

Director 

U.  S.  Naval  Personnel  Research  Activity 

Department  of  the  N i vv 
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Office  of  Naval  Research 
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Commander 

Naval  Ship  Systems  Command  (SHIPS  00V1C) 

Department  of  the  Navy 
Washington,  D.  C.  20360 

Dear  Sir: 

Fig.  7 entitled,  "Three-Dimensional  Projections  of  Two  Quad- 
rature Sampled  Crosscorrelations"  (U)  and  Fig.  10  entitled, 
"Examples  of  Three  Windows  Used  in  Lagged  Cepstrum"  (U)  in 
Quarterly  Progress  Report  No.  7 (U)  under  Contract  NObsr-95181 
were  classified  CONFIDENTIAL  in  error.  Fig.  7 was  numbered 
AS-68-1019  and  Fig.  10  was  numbered  AS-68-1022.  Both  figures 
were  dated  8-1-68.  They  are  UNCLASSIFIED. 


Applied  Research  Laboratories  will  appreciate  your  lining  out 
the  word  "CONFIDENTIAL"  on  these  figures. 
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Technical  Reports  Office 
(Mrs.)  Barbara  M.  Schulze 
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